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The  structures  of  protonated  pentazole  cations  (RNsH+),  oxygen-containing  anions  such  as  N(NC>2)2~,  N03-, 
and  CKV  and  the  corresponding  ion  pairs  are  investigated  by  ab  initio  quantum  chemistry  calculations.  The 
stability  of  the  pentazole  cation  is  explored  by  examining  the  decomposition  pathways  of  several 
monosubstituted  cations  (RN5H+)  to  yield  N2  and  the  corresponding  azidinium  cation.  The  heats  of  formation 
of  these  cations,  which  are  based  on  isodesmic  (bond-type  conserving)  reactions,  are  calculated.  The  proton- 
transfer  reaction  from  the  cation  to  the  anion  is  investigated. 


Introduction 

There  is  growing  interest  in  ionic  liquids  (ILs)  as  solvents 
for  many  applications  including  environmentally  benign  syn¬ 
thesis,  catalysis,  electrochemistry,  and  ion  separation.1-5  The 
chemical  stabilities  and  very  low  vapor  pressures  of  ILs  make 
them  suitable  as  alternatives  to  other  more  volatile  organic 
solvents,  hence  the  name  “green”  solvents.1-3  By  carefully 
selecting  the  starting  cation  and  anion,  the  IL  can  be  customized 
based  on  preselected  characteristics  (e.g.,  moisture  stability, 
viscosity,  density,  miscibility),  thus  allowing  the  needed  flex¬ 
ibility  in  the  design  of  these  materials. 

Despite  the  “designer  solvent”  classification  of  stable  ILs,  it 
is  important  to  remember6  7  that  they  are  not  required  to  be 
stable.  In  fact,  a  substantial  number  of  ILs  are  chemically 
unstable,  releasing  large  amounts  of  energy  in  the  decomposition 
process.  These  so-called  “energetic”  ionic  liquids  may  find  some 
other  useful  applications.8  It  is  interesting  to  note  that  these 
properties  of  ILs  are  analogous  to  the  highly  desirable  charac¬ 
teristics  of  high-energy  density  materials  (HEDMs)  studied 
elsewhere.9-11  The  large  positive  heats  of  formation  of  cyclic 
nitrogen-rich  cations  make  them  promising  precursors  for  the 
synthesis  of  highly  energetic  materials.  One  of  the  most  widely 
studied  IL  cations  is  butyl-methyl-imidazolium  (BMIM+). 
Molecular  dynamics  simulations  of  the  IL  BMIM+/PF6  indicate 
at  least  two  different  time  scales  for  diffusion  in  such  solvents 
and  provide  information  on  preferred  structures.12’13  The  interac¬ 
tion  of  the  surface  layer  of  this  IL  with  water  has  been  studied 
experimentally14  and  the  observed  structural  changes  are 
consistent  with  the  results  of  the  computer  simulations. 

To  increase  the  energy  content  (with  the  side  effect  of 
decreasing  the  stability)  of  the  cation,  BMIM+  may  be  replaced 
by  the  triazolium  and  tetrazolium  cations  with  three  and  four 
nitrogens,  respectively,  in  the  ring.  Recently,  theoretical  studies 
have  been  completed  on  ionic  liquids  based  on  the  triazolium 
(N3C2H4+)  and  tetrazolium  (N4CH3+)  cations,  combined  with 
anions  such  as  dinitramide  (N(N©2)2_),  nitrate  (NO3-),  and 
perchlorate  (C104-).15’16  Another  approach  to  increasing  the 
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(I)  (II)  (III) 

Figure  1.  General  structures  of  monosubstituted  l-H,3-R-l,3-pentazole 
(I),  l-H,2-R-l,2-pentazole  (II),  and  1-H,1-R- 1,1 -pentazole  (III)  cations 
showing  the  atom  numbering  around  the  ring. 

energy  of  the  cations  is  to  introduce  high-energy  substituents 
into  the  ring,  such  as  azido  (— N3)  and  nitro  (— NO2). 

The  present  work  is  a  systematic  investigation  of  the 
structures  and  energetics  of  possible  ionic  liquids  formed  by 
protonated  pentazole  cations,  RNsH+,  and  the  dinitramide, 
nitrate,  and  perchlorate  anions. 

There  have  been  few  attempts  to  synthesize  HNsH+  as  even 
its  parent  precursor,  pentazole  (HN5),  has  proved  to  be  difficult 
to  isolate  and  detect.  The  first  pentazole  compounds  to  be 
synthesized  were  substituted  by  a  phenyl  ring.17-20  Huisgen  and 
Ugi  prepared  their  phenyl- substituted  pentazole  by  reacting  the 
azide  anion,  N3-,  with  diazonium  ions,  ArN2+,  under  careful 
temperature  control.18’19  The  aryl  pentazole  products  readily  lost 
N2,  but  the  resulting  dimethylanilino  pentazole  was  stable 
enough  to  allow  an  X-ray  crystal  structure  determination. 
However,  they  were  not  successful  in  isolating  HN5.  Since  then 
a  significant  goal  of  heterocyclic  chemistry  has  been  to  prepare 
the  parent  pentazole.  In  earlier  attempts,  ozonolytic  degradation 
of  the  aryl  ring  with  a  controlled  amount  of  ozone  was  made 
but  neither  HN5  nor  its  anion  N5-  was  found.21  Recently,  Butler 
et  al.,  working  under  controlled  conditions,  were  able  to  isolate 
and  detect  pentazole  in  solution  as  N5-  with  Zn2+  ion.22 

In  addition  to  the  experimental  investigations,  several  theo¬ 
retical  studies  have  been  performed  on  HN5  as  well.23-27 
Calculations  predict  that  the  half-life  of  HN5  is  only  about  10 
min;21  this  explains  the  difficulty  in  the  direct  detection  of  this 
compound.  Chen23  had  initially  suggested  that  HN5  is  a  strong 
acid,  even  stronger  than  nitric  acid;  hence,  in  solution  a 
significant  portion  is  easily  converted  to  the  anion,  a  fact  finally 
corroborated  by  the  experimental  results  of  Butler  et  al.22  On 
the  other  hand,  attempts  to  isolate  the  protonated  pentazole 
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Figure  2.  Optimized  structures  of  dinitramide,  nitrate,  and  perchlorate 
anions  (lower  structures)  and  their  protonated  counterparts  (nitrogen 
=  blue,  oxygen  =  red,  chlorine  =  green,  carbon  =  black,  hydrogen  = 
white). 

&  3  A. 
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Figure  3.  Monosubstituted  l-H,3-R-pentazole  minimum  energy  struc¬ 
tures  (nitrogen  =  blue,  oxygen  =  red,  fluorine  =  green,  carbon  =  black, 
hydrogen  =  white). 

cation  have  so  far  been  unsuccessful  and  theoretical  calculations 
are  sparse.28’29  The  lack  of  experimental  data  on  protonated 
pentazole  makes  accurate  quantum  chemical  calculations  es¬ 
sential  in  predicting  the  properties  of  these  compounds. 

The  acidity  of  the  azoles  has  been  known  for  quite  sometime 
and  increases  in  the  order  pyrrole  <  pyrazole  <  imidazole  < 
1,2, 3 -triazole  <  1,2,4-triazole  <  tetrazole  <  pentazole.30  This 
means  that  the  protonated  pentazole  cation,  which  is  the  most 
acidic  among  the  protonated  azoles,  is  the  most  difficult  to 
isolate  experimentally.  Mo  et  al.  investigated29  the  basis  set 
effects  on  the  theoretical  (Hartree— Fock)  gas-phase  basicities 
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Figure  4.  Monosubstituted  l-H,2-R-pentazole  cation  minimum  energy 
structures.  Atom  colors  are  defined  in  Figure  3. 

of  several  azoles  including  pentazole  and  its  protonated  form 
with  the  STO-3G,  3-2 1G,  and  6-3 1G  basis  sets,  using  the 
conjugate  gradient  optimization  procedure31  to  obtain  the  starting 
structures.  The  basicity  of  these  compounds  is  predicted  to 
depend  on  the  contributions  of  different  tautomers  of  either  the 
nonprotonated  or  protonated  forms. 

The  monosubstituted  pentazole  cation  has  three  possible 
isomers  (i.e.,  l-H,3-R-l,3-pentazole  (I),  l-H,2-R-l,2-pentazole 
(II),  and  1-H,1-R- 1,1 -pentazole  (III)  ions)  (see  Figure  1).  The 
only  difference  among  the  three  isomers  is  in  the  location  of 
the  substituent  R.  To  find  new  high-energy  IL  precursors  and 
to  obtain  a  more  comprehensive  picture  of  the  stability  of  these 
compounds,  several  topics  will  be  addressed  below.  First,  the 
stability  and  structure  of  several  monosubstituted  pentazole 
cations  are  investigated  by  determining  the  activation  barrier 
and  the  energy  of  decomposition  to  N2  and  the  corresponding 
azidinium  cation.  It  is  expected  that  substituents  may  play  a 
major  role  in  the  thermodynamic  stability  of  the  protonated 
pentazole  ring.  Second,  jr-bonding  and  population  analysis  are 
performed  to  determine  the  amount  of  charge  delocalization 
around  the  ring  in  these  compounds.  One  of  the  characteristics 
of  low  melting  ILs  is  the  presence  of  charge-delocalized  ions. 
Finally,  the  structures  of  dimers  composed  of  the  protonated 
pentazole  and  several  energetic  anions  including  dinitramide, 
nitrate,  and  perchlorate  are  presented.  All  possible  structures, 
including  the  neutral  pairs,  transition  states,  and  ionic  pairs,  are 
described.  The  small  barriers  for  proton  transfer  from  the  cation 
to  the  anion  indicate  that  deprotonation  may  be  a  first  step  in 
the  decomposition  of  these  pentazole-based  ILs. 

Computational  Details 

The  geometries  of  the  monosubstituted  pentazole  cations 
RN5H+  (R  =  H,  F,  CH3,  CN,  NH2,  OH,  OCH3,  CH2N02,  C2H3, 
N3,  NF2),  transition  state  structures  and  decomposition  products 
as  well  as  the  anions  dinitramide,  nitrate,  and  perchlorate  were 
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TABLE  1:  Monosubstituted  l-H,3-R-l,3-Pentazole  (I)  Optimized  Geometries  (Bond  Lengths  A— B  in  A  and  Angles  A— B— C  in 
Degrees)" 


H 

F 

ch3 

CN 

nh2 

OH 

ch2no2 

n3 

nf2 

c2h3 

och3 

N3-R 

1.026 

1.299 

1.476 

1.372 

1.355 

1.322 

1.459 

1.361 

1.463 

1.433 

1.310 

N1-N2 

1.313 

1.318 

1.314 

1.314 

1.322 

1.320 

1.313 

1.320 

1.312 

1.315 

1.321 

N2-N3 

1.313 

1.311 

1.317 

1.325 

1.321 

1.315 

1.317 

1.325 

1.313 

1.323 

1.323 

N3-N4 

1.328 

1.321 

1.334 

1.348 

1.339 

1.331 

1.334 

1.347 

1.332 

1.342 

1.337 

N4-N5 

1.313 

1.319 

1.312 

1.310 

1.314 

1.317 

1.312 

1.314 

1.314 

1.310 

1.320 

N1-N5 

1.328 

1.327 

1.326 

1.328 

1.320 

1.323 

1.327 

1.322 

1.332 

1.326 

1.320 

N1-N2-N3 

99.6 

98.0 

100.5 

99.6 

99.6 

98.9 

100.1 

99.4 

99.2 

100.5 

99.0 

N2-N3-N4 

115.1 

117.4 

113.7 

114.4 

114.5 

115.7 

114.1 

114.4 

115.7 

113.4 

115.2 

N3-N4-N5 

105.1 

103.5 

106.0 

104.9 

105.7 

104.6 

105.6 

105.1 

104.5 

106.0 

104.7 

N4-N5-N1 

105.1 

105.6 

105.0 

105.7 

105.5 

105.3 

105.1 

105.4 

105.4 

105.2 

105.4 

N2-N1-N5 

115.1 

115.4 

114.9 

115.4 

115.2 

115.5 

115.0 

115.7 

115.2 

115.0 

115.7 

N2-N3-R 

122.1 

120.8 

123.3 

121.8 

121.8 

121.7 

122.2 

118.7 

118.0 

122.0 

119.0 

N4-N3-R 

122.8 

121.8 

123.0 

122.8 

122.9 

122.8 

122.8 

126.8 

126.3 

124.6 

125.8 

a  For  numbering  system,  see  Figure  1 . 

TABLE  2:  Monosubstituted  l-H,2-R-l,2-Pentazole  (II)  Optimized  Geometries  (Bond  Lengths  A- 

-B  in  A  and  Angles  A- 

-B-C  in 

Degrees)" 

H 

F 

ch3 

CN 

nh2 

OH 

ch2no2 

n3 

nf2 

c2h3 

och3 

N2-R 

1.027 

1.304 

1.479 

1.369 

1.379 

1.325 

1.474 

1.363 

1.464 

1.434 

1.313 

N1-N2 

1.311 

1.311 

1.315 

1.320 

1.318 

1.312 

1.315 

1.319 

1.311 

1.320 

1.318 

N2-N3 

1.321 

1.309 

1.325 

1.344 

1.330 

1.317 

1.330 

1.334 

1.325 

1.332 

1.322 

N3-N4 

1.327 

1.336 

1.328 

1.322 

1.324 

1.338 

1.326 

1.334 

1.327 

1.326 

1.342 

N4-N5 

1.327 

1.329 

1.324 

1.328 

1.335 

1.321 

1.326 

1.320 

1.334 

1.325 

1.318 

N1-N5 

1.321 

1.327 

1.324 

1.320 

1.314 

1.330 

1.322 

1.329 

1.319 

1.323 

1.331 

N1-N2-N3 

109.3 

112.0 

107.9 

108.7 

108.4 

110.5 

108.5 

109.1 

110.1 

107.6 

110.0 

N2-N3-N4 

105.2 

103.4 

106.2 

104.8 

105.5 

104.4 

105.6 

105.0 

104.4 

106.1 

104.5 

N3-N4-N5 

110.9 

111.5 

110.7 

111.6 

110.9 

111.2 

110.8 

111.3 

111.3 

110.9 

111.3 

N4-N5-N1 

105.2 

105.5 

105.0 

105.4 

105.0 

105.5 

105.3 

105.5 

105.2 

105.0 

105.6 

N2-N1-N5 

109.3 

107.6 

110.2 

109.5 

110.1 

108.4 

109.7 

109.2 

109.0 

110.3 

108.7 

N1-N2-R 

126.3 

122.6 

127.1 

125.5 

121.2 

120.9 

128.1 

120.5 

119.8 

125.8 

120.8 

N3-N2-R 

124.4 

125.3 

125.0 

125.8 

130.4 

128.6 

123.1 

130.3 

130.1 

126.5 

129.2 

"  For  numbering  system,  see  Figure  1 . 


fully  optimized  with  second-order  Mpller-Plesset  perturbation 
theory  (MP2),32’33  using  the  restricted  Hartree— Fock  (RHF)34 
determinant  as  the  reference.  The  completely  renormalized 
coupled-cluster  with  singles,  doubles,  and  noniterative  triples 
(CR— CCSD(T))35-37  method  was  used  to  calculate  single-point 
energies  for  the  cations  at  the  MP2  geometries.  The  basis  set 
used  was  6-311++G(d,p)38’39  for  all  calculations.  Intrinsic 
reaction  coordinate  (IRC)  calculations  were  performed  to 
connect  the  transition  states  with  the  reactants  and  products.40 
A  multiconfiguration  self-consistent  field  (MCSCF)41  calculation 
using  Edmiston-Ruedenberg  molecular  orbitals  (LMO)42  was 
performed  to  determine  the  amount  of  charge  delocalization  in 
the  ring.15’16  Gaussian-2  (G2)43  44  calculations  were  employed 
to  extract  data  for  the  thermochemical  analysis.  Ion  pairs  of 
several  neutral,  transition  state,  and  ionic  dimers  were  optimized 
at  the  RHF  level,  and  MP2  single-point  energies  at  the  RHF 
geometries  were  calculated.  All  results  were  obtained  from  the 
GAMESS45  46  electronic  structure  package,  except  for  the  G2 
results,  which  were  obtained  from  the  Gaussian94  code.47  All 
molecules  were  visualized  with  MacMolPlt.48 

Results  and  Discussion 

Minimum  Energy  Structures.  The  MP2/6-311++G(d,p) 
optimized  structures  for  the  anions  and  their  protonated  coun¬ 
terparts  that  are  used  to  construct  the  dimer  structures  are  shown 
in  Figure  2.  There  is  very  little  change  in  the  geometries  between 
the  protonated  and  deprotonated  structures  of  the  nitrate  and 
perchlorate  species.  However,  the  dinitramide  anion  has  pro¬ 
tonation  sites  at  both  the  central  nitrogen  and  any  of  the  four 
oxygens;  only  the  structure  with  the  preferred  site  of  protonation 
(i.e.,  the  central  nitrogen)  is  shown. 


•*  a 

T  T  -Y 

1  -H,  1  -H~Nj+  +  N2  1  -HJ-F'Nj+  +  N2  FH,l-CHrN3+  +  Nj 


]  -H,  1  -CN-N3+  +  N;  H3N;+  +  2N2  H2N04  +  2N2 


H4CON4  +  2N2 


H4  +  4N2  1-H,l-NF2-N3*  +  Na 


% 


l-Ha-CHjNOj-Nj4  +  N2  H.C^4  +  2N2 


Figure  5.  Isomer  III  minimum  energy  structures  showing  no  stable 
ring  species  from  the  l-H,l-R-pentazole  cation.  See  last  figure  in  Figure 
1  for  the  starting  structures.  Atom  colors  are  defined  in  Figure  3. 
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Figure  6.  Contributing  resonance  structures  for  monosubstituted  1,3- 
pentazole  (IA  and  IB)  and  1,2-pentazole  (IIA  and  IIB). 

The  optimized  geometries  for  the  monosubstituted  pentazole 
cations  I  (Figure  3)  and  II  (Figure  4)  are  summarized  in  Tables 
1  and  2,  respectively.  MP2/6-311++G(d,p)  geometry  optimiza¬ 
tions  of  the  1,1 -pentazole  cations  III  lead  to  the  instantaneous 
decomposition  to  N2  and  the  corresponding  azidinium  cations, 
as  illustrated  in  Figure  5  for  R  =  H— ,  F— ,  CH3— ,  CN— ,  NF2— , 
and  -CH2NO2.  For  R  =  NH2— ,  OH-,  and  -OCH3,  the 
azidinium  cations  break  down  further,  ejecting  another  N2.  For 
R  =  C2H3,  the  1,1 -pentazole  cation  decomposes  into  two  N2 
molecules,  with  cyclization  of  the  remaining  NHC2H3+  species. 
It  is  interesting  to  see  in  Figure  5  that  the  ^-substituted  cation 
spontaneously  decomposes  into  pieces,  forming  N2  molecules 
and  H+. 

Like  their  neutral  parent  HN5,27  the  cationic  pentazole  rings 
in  all  monosubstituted  species  for  isomers  I  and  II  are  planar 
with  almost  equal  N— N  bond  lengths  of  ~1 .32  A  (Tables  1 
and  2).  The  maximum  difference  in  N— N  bond  distances  is 
only  0.038  A  in  isomer  I  for  R  =  CN  with  the  N3— N4  and 
Nl—  N5  bonds  always  slightly  longer  than  the  other  N— N  bonds 
in  the  ring  (see  Figure  1  for  the  atom  numbering).  This 
difference  is  slightly  smaller  in  isomer  II  (about  0.027  A  for  R 
=  F)  where,  in  general,  the  N3— N4  bond  is  a  little  longer  than 
the  remaining  N— N  bonds.  These  results  indicate  that  the 
cationic  pentazole  ring  is  hardly  affected  by  the  attached 
substituents  regardless  of  how  small  (— H)  or  bulky  (— NF2)  the 
groups  are,  and  the  strength  of  their  electron-donating  (—OH 
and  — NH2)  or  electron-withdrawing  (— CN)  character.  The  very 
small  variation  in  the  N— N  distances  also  suggests  significant 
electron  delocalization. 

MP2  geometry  optimization  of  the  nitro  group  (— NO2)  results 
in  a  rather  long  N— N  distance  of  2.684  A,  suggesting  very  weak 
binding  to  the  ring  nitrogen.  This  is  not  surprising,  since  it  was 
shown  previously  that  —  NO2  does  not  bind  well  to  a  nitrogen 
of  the  tetrazolium  ring.13  In  the  present  work,  a  methylene 
spacer,  — CH2NO2,  is  inserted  to  facilitate  binding. 

Relative  Energies  of  Monosubstituted  Pentazole  Cations. 
The  two  stable  ring  isomers  of  the  monosubstituted  pentazole 
cations  (labeled  I  and  II  in  Figure  1)  have  two  possible  resonance 
structures,  shown  in  Figure  6,  that  differ  in  the  positions  of  the 
double  bonds  and  the  location  of  the  cationic  charge.  The 
molecular  and  electronic  structure  of  each  isomer  is  some 
composite  of  these  contributing  resonance  structures  (i.e.,  IA 
and  IB  for  isomer  I  and  IIA  and  IIB  for  isomer  II).  For  all 
substituted  species,  MP2/6-311++G(d,p)  predicts  that  isomer 
I  is  lower  in  energy  than  isomer  II  by  ~16— 22  kcal/mol  (see 
Table  3).  For  comparison,  the  CR— CCSD(T)  single-point 


r 


transition  state 


1-H3-H-N5+ 

Figure  7.  The  MP2/6-311++G(d,p)  relative  energies  in  kcal/mol  for 
proton  rearrangement  from  the  l-H,2-H-pentazole  to  the  1-H,3-H- 
pentazole. 

energies  at  the  MP2  geometries  predict  that  isomer  I  is  lower 
by  ~15 — 21  kcal/mol. 

The  transition  state  for  proton  rearrangement  in  the  hydrogen- 
substituted  cations  (HN5H+)  from  the  1-H,2-H-N5+  species  to 
its  more  stable  1-H,3-H-N5+  counterpart  lies  50  kcal/mol  above 
the  less  stable  isomer  (see  Figure  7).  This  means  that  it  is  rather 
unlikely  for  proton  migration  to  occur,  although  chemical 
intuition  suggests  that  steric  effects  due  to  bulky  substituents 
attached  to  N2  in  the  pentazole  ring,  such  as  —  CH2NO2  and 
— NF2,  may  lower  this  barrier. 

To  assess  the  relative  importance  of  the  contributing  reso¬ 
nance  structures  with  regard  to  the  amount  of  charge  delocal¬ 
ization  around  the  pentazole  ring  in  the  parent  isomers,  an 
MCSCF  jr-LMO  population  and  jr-bond  order  analysis15’16  has 
been  performed.  Recall  that  the  possibility  of  obtaining  a  low 
melting  point  IL  is  enhanced  by  substantial  charge  delocalization 
around  the  protonated  pentazole  ring.  The  diagonal  elements 
of  the  density  matrix  can  be  regarded  as  the  number  of  electrons 
occupying  the  LMO,  while  the  off-diagonal  elements  give  the 
jr-bond  orders.  A  positive  bond  order  reflects  a  bonding 
interaction  and  a  negative  bond  order  indicates  an  antibonding 
interaction.  The  strength  of  these  interactions  is  measured  by 
the  magnitude  of  their  respective  bond  orders. 

The  results  of  the  MCSCF  bonding  analysis  for  the  parent 
pentazole  cations  are  presented  in  Table  4.  For  the  resonance 


TABLE  3:  MP2  and  CR-CCSD(T)/6-311++G(d,p)  Electronic  Energy  Differences,  in  kcal/mol 


H 

F 

ch3 

CN 

nh2 

OH 

ch2no2 

n3 

NF2 

c2h3 

OCR 

MP2 

0.0 

0.0 

0.0 

0.0 

1,3 -Isomer 

0.0  0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CR-CCSD(T) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

MP2 

19.4 

21.7 

18.1 

19.2 

1 ,2-Isomer 

16.7  18.0 

18.0 

15.9 

18.9 

18.0 

16.2 

CR-CCSD(T) 

19.0 

21.1 

17.9 

18.9 

15.3 

17.8 

17.7 

15.7 

18.9 

18.0 

15.9 

Pentazole-Based  Energetic  Ionic  Liquids 
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TABLE  4:  MCSCF/6-311++G(d,p)  it  Orbital  Populations 
and  Bond  Orders 


7i  orbital  populations 


N1 

N2 

N3 

N4 

N5  Total  bonding 

1,3  ion 

1.37 

1.07 

1.39 

1.07 

1.09 

1.89 

1,2  ion 

1.42 

1.42 

1.04 

1.06 

1.04 

1.89 

adjacent  tx  bond  orders 

N1-N2 

N2-N3 

N3-N4 

N4-N5 

N5-N1 

1,3  ion 

0.62 

0.60 

0.52 

0.70 

0.51 

1,2  ion 

0.49 

0.58 

0.62 

0.62 

0.58 

next  neighbor  jx  bond  orders 


N1-N3  N2-N4  N3-N5  N4-N1  N5-N2 


1,3  ion  -0.29  -0.22  -0.17  -0.15  -0.23 

1,2  ion  -0.32  -0.16  -0.04  -0.16  -0.32 

TABLE  5:  Heats  of  Formation  at  298  K  for  Gas-Phase 
Monosubstituted  1,3-  and  1,2-Pentazole  Cations" 


1,3-isomer  1,2-isomer 


R 

AHA 

(IA) 

AH,” 

(IB) 

average 

AH,° 

(IIA) 

AHA 

(IIB) 

average 

H 

285.6 

285.6 

285.6 

4.0 

303.9 

303.9 

303.9 

4.2 

F 

307.9 

307.2 

307.6 

3.4 

328.5 

327.8 

328.2 

3.6 

ch3 

268.3 

267.6 

268.0 

3.1 

284.8 

284.2 

284.5 

3.3 

CN 

358.1 

358.5 

358.3 

3.7 

376.6 

377.0 

376.8 

3.9 

nh2 

302.2 

301.0 

301.6 

3.5 

318.3 

317.1 

317.7 

3.7 

OH 

283.7 

281.8 

282.8 

3.2 

301.2 

299.3 

300.3 

3.4 

CH2N02 

278.8 

277.9 

278.3 

2.1 

296.1 

295.1 

295.6 

2.2 

n3 

393.6 

391.0 

392.3 

3.5 

409.2 

406.6 

407.9 

3.6 

nf2 

317.3 

313.9 

315.6 

2.6 

335.4 

332.1 

333.7 

2.7 

c2h3 

294.6 

293.3 

294.0 

3.0 

311.1 

309.8 

310.5 

3.2 

och3 

278.1 

276.4 

277.2 

2.7 

293.9 

292.2 

293.0 

2.9 

a  IA,  IB,  IIA,  and  IIB  refer  to  particular  resonance  hybrids  with  their 
average  values  being  the  predicted  heats  (the  average  is  given  in  both 
kcal/mol  and  kcal/g). 


structures  of  isomer  I,  if  IA  dominated  over  IB,  N1  would  have 
two  jr-electrons  and  all  other  atoms  in  the  ring  would  possess 
one  jr-electron  (see  Figure  6).  Table  4  shows  that  N1  and  N3 
have  jr-populations  of  1.37  and  1.39  electrons,  respectively, 
indicating  nearly  equal  contributions  from  the  two  resonance 
structures.  So,  these  two  nitrogens  equally  share  the  formal 
charge  of  +1  in  isomers  IA  and  IB.  For  isomer  II,  both  N1  and 
N2  have  the  same  population  of  1 .42  jr-electrons,  which  means 
that  the  ring  is  a  hybrid  between  IIA  and  IIB  and  that  the 
resonance  structures  are  of  equal  importance.  These  observations 
are  further  validated  by  the  calculated  jr-bond  orders  of  0.50— 
0.70  between  adjacent  atoms.  These  bond  orders  between  0  and 
1  also  suggest  that  the  ring  is  delocalized.  For  both  isomers,  all 
next  nearest  neighbor  interactions  are  antibonding.  The  calcu¬ 
lated  total  jr-bonding  of  1.89,  obtained  by  summing  all 
jr-bonding  and  all  jr-antibonding  interactions,  is  very  close  to 
the  total  bonding  of  2  jr-bonds  in  diagrams  I  A,  IB,  IIA,  and 
IIB  shown  in  Figure  6. 

Heats  of  Formation.  The  amount  of  energy  contained  in 
the  monosubstituted  pentazole  cations  may  be  quantified 
via  their  heats  of  formation,  calculated  here  by  combining 
their  MP2/6-311++G(d,p)  electronic  energies  with  the  data 
obtained  from  the  G2  calculations.  It  is  expected  that  the 
pentazole  cation  in  which  all  the  atoms  in  the  ring  are  nitrogens 
will  have  a  larger  heat  of  formation  than  tetrazolium  or 
triazolium. 

The  calculated  gas-phase  heats  of  formation  for  both 
monosubstituted  pentazole  cations  are  presented  in  Table  5. 


u  A 

4  A 

Figure  8.  Monosubstituted  l-H,3-R-pentazole  transition  states.  Atom 
colors  are  defined  in  Figure  3. 

To  obtain  accurate  heats  of  formation  for  these  cations,  sever¬ 
al  isodesmic  (bond-type  conserving)  reactions  are  utilized. 
Because  there  are  two  contributing  resonance  structures  each 
for  both  isomers,  four  isodesmic  reactions  are  considered  (two 
for  isomer  I,  i.e.,  IA  and  IB,  and  two  for  isomer  II,  i.e.,  IIA  and 
IIB): 

IA  +  7  NH3  —  2  H2N-NH2  +  2  HN=NH  +  H2N— R  +  H4N+ 
IB  +  7  NH3  —  2  H2N-NH2  +  2  HN=NH  +  H3N  +  H3N-R+ 
IIA  +  7  NH3  —  2  H2N-NH2  +  2  HN=NH  +  H2N— R  +  H4N+ 
IIB  +  7  NH3  —  2  H2N-NH2  +  2  HN=NH  +  H3N  +  H3N-R+ 

Because  isodesmic  reactions  conserve  the  number  of  formal 
bond  types  on  both  sides  of  a  reaction,  the  correlation  energy49 
error  is  minimized,  and  hence  the  computed  enthalpies  of 
formation  are  more  accurate  because  of  error  cancellations. 
Table  5  shows  that  the  calculated  heats  of  formation  for  IA  and 
IB  (or  IIA  and  IIB)  are  almost  always  equal  with  the  largest 
deviation  of  only  3.3  kcal/mol  obtained  from  the  NF2-substituted 
ion.  Thus,  the  average  values  of  IA  and  IB  and  of  IIA  and  IIB 
can  be  regarded  as  the  heats  of  formation  for  the  monosubsti¬ 
tuted  pentazole  cations  I  and  II,  respectively. 

As  expected,  the  average  heats  of  formation  for  the  mono¬ 
substituted  pentazole  cations  are  large,  about  4.1  kcal/g  for  the 
parent  H-substituted  pentazole  cations.  In  comparison,  the 
predicted  heats  of  formation  for  the  parent  triazolium  and 
tetrazolium  cations  are  2.8 15  and  3.6 16  kcal/g,  respectively.  This 
is  further  evidence  that  these  pentazole  ions  are  relatively 
unstable  and,  therefore,  will  likely  be  difficult  to  synthesize, 
isolate,  and  detect  much  like  their  parent  neutral  pentazole.  The 
heats  of  formation  of  the  monosubstituted  pentazole  cations  I 
are  always  smaller  by  about  15.5—20.6  kcal/mol  relative  to 
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TABLE  6:  Monosubstituted  l-H,3-R-l,3-Pentazole  Transition  State  Geometries  (Bond  Lengths  A— B  in  A  and  Angles  A— B— C 
in  Degrees)" 


H 

F 

ch3 

CN 

nh2 

OH 

CH2NO2 

n3 

nf2 

c2h3 

och3 

N3-R 

1.028 

1.303 

1.468 

1.349 

1.305 

1.310 

1.449 

1.330 

1.430 

1.408 

1.288 

N1-N2 

1.243 

1.264 

1.259 

1.251 

1.288 

1.274 

1.246 

1.289 

1.247 

1.272 

1.282 

N2-N3 

1.243 

1.267 

1.238 

1.253 

1.275 

1.267 

1.249 

1.270 

1.259 

1.246 

1.272 

N3-N4 

1.699 

1.778 

1.761 

1.793 

1.789 

1.786 

1.704 

1.815 

1.739 

1.815 

1.793 

N4-N5 

1.181 

1.181 

1.182 

1.180 

1.193 

1.186 

1.182 

1.191 

1.184 

1.183 

1.189 

N1-N5 

1.699 

1.586 

1.643 

1.602 

1.523 

1.558 

1.684 

1.511 

1.634 

1.584 

1.546 

N1-N2-N3 

115.8 

112.0 

116.3 

114.7 

111.4 

112.2 

115.7 

111.7 

113.7 

115.5 

111.8 

N2-N3-N4 

106.3 

105.5 

103.8 

103.6 

103.7 

104.3 

105.4 

103.2 

105.6 

102.4 

104.2 

N3-N4-N5 

105.7 

100.9 

104.4 

102.6 

101.2 

101.4 

106.2 

100.4 

103.6 

102.3 

101.1 

N4-N5-N1 

105.7 

110.2 

107.4 

109.0 

110.5 

110.0 

105.5 

111.3 

107.6 

109.7 

110.5 

N2-N1-N5 

106.3 

110.4 

107.6 

109.9 

111.3 

110.9 

106.9 

111.9 

109.1 

109.3 

111.1 

N2-N3-R 

125.0 

120.5 

129.8 

129.0 

124.9 

121.0 

126.0 

122.0 

120.2 

129.3 

122.9 

N4-N3-R 

125.4 

126.8 

124.2 

125.8 

128.6 

127.6 

120.9 

129.2 

129.4 

126.2 

127.4 

a  For  numbering  system,  see  Figure  1 . 


cations  II.  These  are  all  larger  than  the  inferred  heats  of 
formation  of  the  oxadiaziridine-based  HEDMs.11 

Of  course,  the  nature  of  the  substituent  attached  to  the 
pentazole  ring  affects  the  computed  enthalpy  of  formation.  By 
replacing  the  —  H  substituent  in  the  parent  pentazole  cation  by 
electron-donating  groups  (EDGs),  the  calculated  heats  of 
formation  decrease,  while  electron-withdrawing  groups  (EWGs) 
tend  to  increase  the  enthalpies.  For  the  monosubstituted  isomer 
I,  the  heats  of  formation  (in  kcal/mol)  are  decreased  by  the 
presence  of  EDGs  by  as  little  as  2.8  (—OH)  to  as  much  as  17.6 
(— CH3)  and  increased  by  EWGs  by  as  little  as  22.0  (— F)  to  as 
much  as  106.8  (— N3)  relative  to  the  enthalpy  for  the  parent 
H-substituted  1,3-pentazole  ion  (285.6  kcal/mol).  These  differ¬ 
ences  in  enthalpies  are  3.7,  19.4,  24.2,  and  104.0  kcal/mol  for 
the  —OH,  — CH3,  — F,  and  — N3  substituents,  respectively, 
relative  to  the  parent  1,2-pentazole  ion  (303.9  kcal/mol). 
However,  the  cations  with  nitrogen- containing  EDGs  such  as 
— NH2  and  —  NF2  have  higher  heats  of  formation  in  opposition 
to  the  main  trend.  In  these  instances,  the  resulting  increase  in 
the  calculated  enthalpy  is  due  to  the  additional  nitrogen  adjacent 
to  the  pentazole  ring.  All  listed  heats  of  formation  in  Table  5 
are  large  and  positive,  but  the  species  with  the  largest  enthalpies 
on  a  per  weight  basis  include  the  parent  pentazole  cations  with 
heats  of  formation  of  about  4.1  kcal/g  and  the  CN-substituted 
ions  with  energies  of  about  3.8  kcal/g. 

Transition  States  and  Decomposition  Products.  The  transi¬ 
tion  state  structures  for  the  decomposition  of  the  monosubstituted 
cations  I,  RNsH+  — ►  RN3H+  +  N2,  (Figure  8)  show  a 
lengthening  of  the  Nl—  N5  and  N3— N4  bonds  by  about  0.2— 
0.5  A  and  shortening  of  the  N4— N5  distance  by  about  0.13  A 
(see  Figure  1  and  Table  6),  in  the  process  of  eliminating  an  N2 
molecule.  In  general,  the  N3— N4  bond  length  is  longer  than 
the  Nl— N5  distance,  implying  that  the  transition  state  is 
obtained  via  an  asymmetric  ring  opening.  The  exception  is  the 
parent  1,3-pentazole  ion  (i.e.,  the  H-substituted)  formed  through 
a  symmetric  ring  opening  mechanism  with  equal  N3— N4  and 
Nl—  N5  bond  lengths.  On  the  other  hand,  the  decomposition  of 
the  monosubstituted  cations  II  involves  two  steps.  First,  the  N2— 
N3  distance  increases  by  ~0.5— 0.7  A,  while  the  N3— N4  bond 
decreases  by  M).17  A  forming  ring-opened  transition  states  (see 
Figure  9  and  Table  7),  leading  to  the  local  minimum  structures 
shown  in  Figure  10  and  quantified  in  Table  8.  There  is  a  slight 
increase  in  the  N4— N5  bond  length  of  about  0.05  A,  which 
indicates  that  a  succeeding  step  to  form  N2  is  imminent.  Indeed, 
the  second  step  involves  the  formation  of  transition  states  (Figure 
1 1)  in  which  the  N4— N5  distance  is  increased  further  by  about 
0.6  A  (Table  9).  It  is  interesting  to  note  from  the  structures  in 
Figure  1 1  that  the  hydrogen  in  the  azidinium  cation  located  in 


the  Nl  position  is  oriented  in  such  a  way  that  it  may  be  involved 
in  a  rearrangement  reaction  into  the  more  stable  form  in  which 
the  hydrogen  is  farthest  away  from  the  substituents.  The  one 
exception  to  this  is  the  CH2NO2- substituted  cation  in  which  the 
hydrogen  in  the  N 1  position  is  transferred  to  one  of  the  oxygens 
in  the  CH2NO2  moiety. 

The  structures  of  the  decomposition  products  of  the  various 
monosubstituted  pentazole  cations  for  both  isomers  are  shown 
in  Figures  12  and  13.  The  optimized  geometries  for  the  products 
are  listed  in  Tables  10  and  11.  The  decomposition  of  the 
monosubstituted  pentazole  cations  involves  the  formation  of  N2. 
The  resulting  bond  distance  of  1.12  A  in  the  product  N— N  bonds 
for  all  isomer  I  and  II  decomposition  products  is  in  good 
agreement  with  the  experimental  value  of  1.098  A.  Figure  13 
confirms  that  the  hydrogen  in  the  azidinium  cation,  originally 
located  in  the  Nl  position,  migrates  to  form  the  more  stable 
isomer,  similar  to  that  of  the  isomer  I  decomposition  products. 
As  noted  above,  the  exception  is  the  CH2N02-substituted 
azidinium  ion  in  which  the  hydrogen  is  transferred  to  the  oxygen 


f 


*7 

Figure  9.  Monosubstituted  l-H,2-R-pentazole  ring-opened  transition 
states.  Atom  colors  are  defined  in  Figure  3. 
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TABLE  7:  Monosubstituted  l-H,2-R-l,2-Pentazole  Ring-Opened  Transition  State  Geometries  (Bond  Lengths  A— B  in  A  and 
Angles  A— B— C  in  Degrees)" 


H 

F 

ch3 

CN 

nh2 

OH 

CH2NO2 

n3 

nf2 

c2h3 

och3 

N2-R 

1.029 

1.330 

1.470 

1.347 

1.353 

1.345 

1.460 

1.364 

1.461 

1.402 

1.327 

N1-N2 

1.285 

1.318 

1.275 

1.301 

1.310 

1.316 

1.288 

1.313 

1.303 

1.284 

1.319 

N2-N3 

2.012 

1.836 

2.065 

1.985 

2.013 

1.929 

2.033 

1.969 

1.956 

2.029 

1.948 

N3-N4 

1.152 

1.171 

1.149 

1.156 

1.156 

1.163 

1.152 

1.160 

1.157 

1.153 

1.162 

N4-N5 

1.368 

1.390 

1.355 

1.373 

1.365 

1.377 

1.375 

1.365 

1.387 

1.352 

1.369 

N1-N5 

1.318 

1.309 

1.330 

1.318 

1.320 

1.313 

1.315 

1.319 

1.307 

1.332 

1.317 

N1-N2-N3 

93.1 

96.7 

92.0 

93.6 

92.7 

94.2 

92.8 

93.4 

93.9 

92.6 

93.9 

N2-N3-N4 

89.6 

95.0 

88.0 

90.9 

89.9 

92.6 

89.1 

91.4 

91.9 

89.2 

91.8 

N3-N4-N5 

132.5 

126.5 

134.5 

131.1 

132.3 

129.4 

132.8 

130.8 

129.9 

133.4 

130.5 

N4-N5-N1 

101.3 

101.9 

101.0 

101.9 

101.6 

101.7 

101.4 

101.7 

101.5 

101.3 

101.6 

N2-N1-N5 

123.5 

119.7 

124.4 

122.4 

123.0 

121.9 

123.8 

122.3 

122.8 

123.5 

121.9 

N1-N2-R 

113.6 

113.8 

119.1 

120.3 

115.6 

111.6 

118.2 

112.6 

109.6 

120.7 

112.7 

N3-N2-R 

153.3 

147.7 

149.0 

146.1 

150.1 

148.8 

148.7 

149.7 

156.5 

146.7 

149.5 

a  For  numbering  system,  see  Figure  1 . 


in  the  nitro  group  instead,  producing  three  decomposition 
products:  N2,  HONO,  and  H2CN3+. 

The  activation  energies  for  the  decomposition  reactions  are 
listed  in  Table  12.  Isomer  I  decomposes  in  a  single  step  via  an 
asymmetrically  ring-opened  transition  state  (see  Figure  8  and 
Table  6),  while  isomer  II  decomposes  in  two  steps  via  a  ring 
opening  followed  by  a  subsequent  proton  migration  (Figures  9 
and  11  and  Tables  7  and  9).  All  of  the  monosubstituted  isomer 
I  species  have  CR— CCSD(T)  activation  energies  in  the  20—30 
kcal/mol  range.  The  CR— CCSD(T)  £ai  activation  energies  of 
isomer  II  are  larger  by  ~2— 7  kcal/mol,  except  for  R  =  F,  CN, 
and  NF2.  Previous  studies  on  several  neutral  azoles,  including 
triazole,  tetrazole,  and  pentazole,  show  that  the  activation 
energies  for  this  family  of  azoles  depend  on  the  number  of  ring 
nitrogen  atoms  and  increase  in  the  order  pentazole  <  tetrazole 
<  triazole.27  It  is  expected  that  their  monosubstituted  cationic 
counterparts  will  follow  the  same  trend. 


Figure  10.  Local  minima  structures  connecting  the  transition  states 
of  the  monosubstituted  l-H,2-R-pentazole  cations.  Atom  colors  are 
defined  in  Figure  3. 


Table  12  shows  that  EDGs  in  isomer  I  such  as  —  CH3,  —  CH2- 
N02,  and  —  C2H3  tend  to  give  the  largest  activation  energies 
(24.0,  22.6.  and  21.8  kcal/mol,  respectively)  while  EWGs  such 
as  — N3  (16.0  kcal/mol)  and  — CN  (19.0  kcal/mol)  give  the 
smallest  barriers.  The  £a  1  activation  energies  are  slightly  larger 
in  isomer  II,  but  the  observed  trend  is  the  same.  Isomer  II  with 
EDGs  attached  (e.g.,  -CH3  (37.0  kcal/mol)  and  -C2H3  (32.0 
kcal/mol))  possess  larger  activation  energies  than  their  EWG- 
substituted  counterparts  such  as  — N3  (23.4  kcal/mol)  and  — CN 
(19.0  kcal/mol).  In  general,  the  CR— CCSD(T)  activation  barriers 
are  at  most  5  kcal/mol  larger  than  the  MP2  activation  energies. 
The  inductive  effects  (e.g.,  in  —  CH3  and  —  CH2N02)  and 
resonance  donating  effects  (e.g.,  in  —  C2H3)  brought  about  by 
the  EDGs  help  to  stabilize  the  ring  by  increasing  the  electron 
density  in  the  ring,  thereby  resulting  in  higher  activation 
energies.  Although  isomer  I  is  lower  in  energy,  isomer  II  is 
more  stable  to  decomposition. 

An  important  factor  in  the  design  of  high-energy  materials 
is  the  energy  released  on  decomposition.  The  performance  of  a 
highly  energetic  system  is  greatly  enhanced  with  an  increase  in 
the  released  energy.  The  decomposition  reactions  for  all 
monosubstituted  species  are  exothermic  (see  Table  13).  In 
general,  the  decomposition  of  isomer  II  is  more  exothermic  by 
at  least  14  kcal/mol  relative  to  isomer  I.  As  for  the  activation 
energy,  the  decomposition  energy  is  greatly  affected  by  the 
substituent  attached  in  the  cationic  pentazole  rings  of  the  two 
isomers.  The  decomposition  becomes  more  exothermic  with  an 
increase  in  the  electron-withdrawing  character  of  the  sub¬ 
stituent.  The  highest  decomposition  energies  are  found  in  the 
NH2-  and  N3- substituted  cations  for  isomer  I  (—24.5  and 
—24.0  kcal/mol,  respectively)  and  the  N3-  and  CH2N02- 
substituted  cations  for  isomer  II  (—58.9  and  —56.1  kcal/mol, 
respectively).  The  observed  trend  in  which  a  more  exo¬ 
thermic  reaction  is  accompanied  by  a  lower  activation  barrier 
seems  to  follow  Hammond’s  postulate,50  according  to  which  a 
lower  activation  energy  corresponds  to  a  more  exothermic 
reaction.  By  combining  the  energy  profiles  for  the  decompo¬ 
sition  of  the  monosubstituted  pentazole  cations  including  the 
minimum  energies,  the  activation  energies,  and  the  energies 
for  decomposition,  a  general  potential  energy  pathway  for 
decomposition  can  be  deduced  for  isomers  I  and  II  (Figure 
14).  Finally,  note  that  there  are  several  cases  in  both  Tables 
12  and  13  for  which  there  is  significant  (>10  kcal/mol)  dis¬ 
agreement  between  the  MP2  and  CR— CCSD(T)  energy  differ¬ 
ences.  For  example,  the  two  levels  of  theory  differ  by  ~14 
kcal/mol  in  their  prediction  of  the  net  reaction  energy  for 
the  parent  species  (Table  13),  with  R  =  H.  Such  disagreement 
is  often  indicative  of  some  nontrivial  configurational  mixin 
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TABLE  8:  Monosubstituted  l-H,2-R-l,2-Pentazole  Geometries  for  the  Intermediates  Connecting  the  Two  Transition  States 
(Bond  Lengths  A— B  in  A  and  Angles  A— B— C  in  Degrees)" 


H 

F 

ch3 

CN 

nh2 

OH 

CH2N02 

n3 

nf2 

c2h3 

och3 

N2- 

-R 

1.029 

1.345 

1.460 

1.353 

1.299 

1.334 

1.462 

1.344 

1.466 

1.389 

1.315 

Nl- 

-N2 

1.265 

1.287 

1.254 

1.290 

1.270 

1.275 

1.269 

1.279 

1.281 

1.268 

1.280 

N2- 

-N3 

3.295 

3.353 

3.306 

3.305 

3.297 

3.346 

3.302 

3.323 

3.338 

3.306 

3.328 

N3- 

-N4 

1.135 

1.137 

1.136 

1.137 

1.139 

1.138 

1.137 

1.139 

1.136 

1.137 

1.139 

N4- 

-N5 

1.307 

1.305 

1.299 

1.309 

1.281 

1.293 

1.305 

1.288 

1.312 

1.294 

1.289 

Nl- 

-N5 

1.329 

1.329 

1.342 

1.320 

1.377 

1.346 

1.329 

1.351 

1.318 

1.345 

1.353 

Nl- 

-N2- 

-N3 

75.5 

74.2 

75.1 

75.5 

76.2 

74.4 

75.6 

75.1 

74.2 

75.3 

74.9 

N2- 

-N3- 

-N4 

38.0 

37.6 

38.7 

38.3 

40.3 

38.3 

38.4 

39.1 

37.6 

39.3 

39.1 

N3- 

-N4- 

-N5 

169.3 

168.2 

169.9 

168.7 

170.8 

169.1 

169.2 

169.8 

168.2 

170.0 

169.7 

N4- 

-N5- 

-Nl 

109.7 

110.7 

109.6 

110.4 

110.3 

110.9 

110.3 

110.8 

110.3 

110.1 

110.7 

N2- 

-Nl- 

-N5 

125.4 

125.4 

126.1 

124.6 

124.0 

125.4 

124.9 

124.8 

126.1 

125.4 

125.0 

Nl- 

-N2- 

-R 

108.1 

106.6 

114.6 

114.0 

118.2 

108.4 

116.3 

110.2 

103.7 

115.8 

109.5 

N3- 

-N2- 

-R 

176.4 

179.1 

170.3 

170.5 

165.6 

177.2 

168.0 

174.8 

178.0 

169.0 

175.5 

a  For  numbering  system,  see  Figure  1 . 


g  in  which  case  the  coupled  cluster  method  is  generally  more 
reliable. 

Ion  Pair  Interactions.  The  RHF/6-311++G(d,p)  optimized 
structures  of  the  neutral  and  ionic  pairs,  as  well  as  the  transition 
states  that  connect  them,  for  all  dimers  formed  by  combining 
either  of  the  H- substituted  cationic  pentazole  isomers  with  the 
three  selected  anions  are  shown  in  Figures  15—23.  All  MP2/ 
6-311++G(d,p)  optimizations  and  structure  determinations  were 
unsuccessful  in  locating  ionic  dimers;  only  hydrogen-bonded 
neutral  pairs  seem  to  exist.  However,  MP2  single-point  energies 
at  the  RHF  geometries,  referred  to  as  MP2//RHF,  are  calculated 
for  comparison. 

The  results  of  the  RHF/6-311++G(d,p)  geometry  optimiza¬ 
tions  on  several  different  orientations  of  either  of  the  two 
minimum  energy  protonated  pentazole  isomers  and  the  mini¬ 
mum  energy  dinitramide  anion  indicate  that  the  dimer  structure 
in  Figure  15,  brought  about  by  the  transfer  of  the  proton  from 
the  1,2-isomer  to  the  dinitramide,  is  the  most  stable  neutral 
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Figure  11.  Monosubstituted  l-H,2-R-pentazole  structures  for  the 
second  transition  state.  Atom  colors  are  defined  in  Figure  3. 


dimer.  Henceforth,  this  structure  is  taken  as  the  zero  energy 
structure  against  which  all  other  energies  of  the  remaining 
protonated  pentazole— dinitramide  dimer  structures  are  com¬ 
pared.  However,  the  neutral  dimer  structure  formed  by  the 
protonated  1,3-pentazole  and  dinitramide  is  higher  in  energy 
by  only  0.4  and  1.0  kcal/mol  at  the  RHF  and  MP2  levels  of 
theory,  respectively.  The  MP2//RHF  energy  is  taken  to  be  the 
MP2  zero  of  energy. 

In  the  search  for  the  minimum  energy  protonated  pentazole— 
dinitramide  dimer  structures,  several  RHF  optimizations  were 
performed.  The  starting  structures  were  varied  with  an  adequate 
cation— anion  separation  to  remove  any  biases  toward  the  final 
optimized  structures.  Figures  16,  18,  19,  and  21  show  that  only 
a  few  neutral  and  ionic  dimer  pairs  out  of  about  20  starting 
structures  were  found,  depending  on  the  location  of  the 
protonation  in  dinitramide.  In  particular,  the  number  of  unique 
protonated  pentazole— dinitramide  ionic  dimer  structures  is  only 
four  for  isomer  I  (Figure  16)  and  three  for  isomer  II  (Figure 
19).  The  small  number  of  dimer  structures  found  is  due  to  the 
highly  symmetric  character  of  the  starting  protonated  pentazole 
and  that  the  cationic  pentazole  ring  is  composed  of  all  nitrogens, 
which  means  that  some  of  the  initial  dimer  orientations  converge 
to  the  same  final  structure. 

The  corresponding  neutral  dimer  structures  resulting  from 
proton-transfer  reactions  from  the  cation  to  the  anion  are  shown 
in  Figures  18  and  21.  The  transition  states  linking  the  ionic  pairs 
to  the  neutral  dimers  are  shown  in  Figures  17  and  20.  As  noted 
previously  (Figure  16),  the  common  zero  of  energy  for  all 
structures  is  the  most  stable  pentazole— dinitramide  neutral 
dimer,  shown  here  as  the  leftmost  structure  in  Figure  21.  The 
calculated  MP2  activation  energies  show  shallow  barriers  for 
proton  transfer  with  the  highest  value  only  ~5.6  kcal/mol.  This 
is  indicative  of  the  very  short  lifetimes  of  these  ionic  pairs  if 
they  were  synthesized.  In  fact,  only  two  of  the  ionic  dimers 
have  a  barrier;  in  almost  all  cases,  the  proton  is  spontaneously 
transferred  to  the  dinitramide  anion  to  form  the  neutral  pairs. 

The  calculated  energies  for  the  neutral  dimer  structures  in 
Figures  18  and  21  show  that  the  preferred  protonation  site  in 
the  dinitramide  anion  is  indeed  the  central  nitrogen.  For  all  seven 
neutral  structures  of  the  pentazole— dinitramide  pairs,  the  energy 
separation  is  small  with  a  maximum  difference  at  the  RHF  and 
MP2  levels  of  only  11.9  and  14.0  kcal/mol,  respectively.  The 
structures  in  Figures  18  and  21  show  that  the  predominant  form 
of  interaction  is  intermolecular  hydrogen  bonding,  either  the 
N— H— N  or  O— H— N  type,  but  a  few  dimer  structures  possess 
intramolecular  hydrogen  bonding  within  the  dinitramic  acid 
moiety  such  as  the  O— H— O  bonding  interaction  found  in  Figure 
18(d). 
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TABLE  9:  Monosubstituted  l-H,2-R-l,2-Pentazole  Geometries  for  the  Second  Transition  State  (Bond  Lengths  A— B  in  A  and 
Angles  A— B— C  in  Degrees)" 


H 

F 

ch3 

CN 

nh2 

OH 

CH2N02 

n3 

nf2 

c2h3 

och3 

N2-R 

1.029 

1.359 

1.482 

1.356 

1.303 

1.374 

1.439 

1.220 

1.490 

1.422 

1.296 

N1-N2 

1.255 

1.308 

1.240 

1.270 

1.249 

1.312 

1.291 

1.294 

1.272 

1.248 

1.273 

N2-N3 

3.925 

4.123 

3.945 

3.957 

2.291 

3.844 

3.469 

4.078 

4.151 

3.973 

4.326 

N3-N4 

1.122 

1.122 

1.122 

1.122 

1.144 

1.120 

1.130 

1.140 

1.122 

1.122 

1.140 

N4-N5 

1.961 

2.113 

1.939 

2.000 

1.392 

2.583 

1.524 

1.290 

2.095 

1.965 

1.289 

N1-N5 

1.217 

1.214 

1.227 

1.213 

1.400 

1.217 

1.253 

1.402 

1.206 

1.226 

1.378 

N1-N2-N3 

67.2 

65.4 

65.7 

67.3 

97.6 

77.9 

75.8 

43.2 

62.8 

65.5 

27.4 

N2-N3-N4 

33.8 

25.0 

34.2 

33.4 

43.9 

33.5 

33.6 

34.3 

32.8 

34.4 

10.4 

N3-N4-N5 

161.0 

162.6 

160.1 

161.7 

134.1 

160.7 

154.2 

167.7 

164.1 

160.5 

167.9 

N4-N5-N1 

102.1 

102.6 

102.1 

102.1 

110.0 

81.1 

109.0 

108.8 

102.7 

102.1 

109.5 

N2-N1-N5 

146.8 

149.7 

148.1 

147.1 

113.2 

155.7 

123.9 

116.3 

153.6 

148.5 

116.2 

N1-N2-R 

111.2 

107.7 

118.3 

118.3 

119.9 

109.8 

111.0 

120.9 

104.4 

118.7 

110.0 

N3-N2-R 

178.3 

165.8 

175.9 

174.4 

142.0 

109.4 

170.8 

149.3 

167.2 

175.0 

131.8 

a  For  numbering  system,  see  Figure  1 . 


The  interactions  between  each  of  the  protonated  pentazole 
isomers  with  the  nitrate  and  the  perchlorate  anions  also  were 
examined.  The  MP2  optimizations  result  in  an  instantaneous 
proton  transfer  from  the  protonated  pentazole  to  one  of  the 
oxygens  of  the  anion,  whereas  RHF  optimizations  predict  that 
both  neutral  and  ionic  dimers  are  local  minima.  So,  electron 
correlation  is  clearly  important  in  determining  which  structures 
exist.  Figures  22  and  23  show  the  neutral  and  ionic  dimer 
structures  and  the  transition  states  that  connect  them.  The  pair 
comprising  the  H-substituted  1,2-pentazole  and  either  the  nitrate 
or  perchlorate  anion  is  predicted  to  be  lower  in  energy  than  the 
interaction  between  the  H-substituted  1,3-pentazole  and  these 
anions.  The  zero  energy  structure  for  the  neutral  pentazole— 
nitrate  pair,  resulting  from  the  transfer  of  a  proton  from  isomer 
II  to  the  nitrate,  is  shown  in  the  bottom  right  corner  of  Figure 
22.  The  lowest  energy  protonated  1,3-pentazole— nitrate  pair 
obtained  in  the  MP2  (RHF)  calculation  is  only  1.7  kcal/mol 
(1.3  kcal/mol)  higher  in  energy  relative  to  this  structure.  In  the 
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protonated  pentazole— perchlorate  pair,  the  MP2  energy  separa¬ 
tion  between  the  two  neutral  dimer  structures  is  2.1  kcal/mol. 

Finally,  the  MP2  heats  of  reaction  for  the  formation  of  either 
the  neutral  or  ionic  dimers  from  the  individual  ions  are  included 
as  the  third  number  in  Figures  16,  18,  19,  21,  22,  and  23.  The 
mixing  of  the  protonated  pentazole  with  any  of  the  three  anions 
produces  either  neutral  or  ionic  pairs,  although  MP2  optimiza¬ 
tions  result  in  proton  transfer.  These  reactions  are  all  highly 
exothermic.  As  expected,  the  heats  of  reaction  for  the  formation 
of  the  neutral  pairs  are  always  larger  (more  exothermic)  than 
those  for  the  ionic  pairs,  by  at  least  16  kcal/mol  in  the  protonated 
pentazole— dinitramide,  31  kcal/mol  in  the  protonated  penta¬ 
zole— nitrate,  and  20  kcal/mol  in  the  protonated  pentazole— 
perchlorate  systems.  The  neutral  dimer  containing  1,3-pentazole 


is  lower  in  energy  than  the  corresponding  1,2-isomer  by  only 
1.2,  1.7,  2.2  kcal/mol  for  the  dinitramide,  nitrate,  and  perchlorate 

: 

1-R3-H-N/  +  N; 

1-H.3-F-NS*  +  N, 

t 

l-H.3-CHrNr  +  N: 

: 

l-H53-CN-N3+  +  N2 

X 

1  -H,3-NHrN/  +  N, 

: 

FH3-0H-V  +  N2 

: 

1-H,3-OCH3Nj+  +  N2 

: 

jTS 

HIV  +  N2 

: 

1-H3-NF2-N3+  +  N2 

9*  /  ' 

1  -H,3-CH;NOj-N3*  +  N,  ]  -H.3-C,Hj-N}*  +  N, 

Figure  12.  Monosubstituted  l-H,3-R-pentazole  decomposition  prod¬ 
ucts.  Atom  colors  are  defined  in  Figure  3. 


: 

HjCN^  +  HONO  +  N* 

Figure  13.  Monosubstituted  l-H,2-R-pentazole  decomposition  prod¬ 
ucts.  Atom  colors  are  defined  in  Figure  3. 


t 

1-H,3-C2Hj-Nj+~h  N2 
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TABLE  10:  Monosubstituted  l-H,3-R-l,3-Pentazole  Decomposition  Product  Geometries  (Bond  Lengths  A— B  in  A  and  Angles 
A— B— C  in  Degrees)" 


H 

F 

ch3 

CN 

nh2 

OH 

CH2NO2 

n3 

nf2 

c2h3 

och3 

N3-R 

1.026 

1.348 

1.479 

1.336 

1.247 

1.346 

1.467 

1.362 

1.522 

1.429 

1.329 

N1-N2 

1.185 

1.179 

1.194 

1.189 

1.247 

1.196 

1.183 

1.207 

1.191 

1.201 

1.203 

N2-N3 

1.185 

1.226 

1.182 

1.192 

1.379 

1.205 

1.174 

1.198 

1.196 

1.162 

1.210 

N3-N4 

3.000 

3.052 

2.920 

2.943 

3.899 

2.977 

3.221 

2.968 

3.953 

3.094 

2.968 

N4-N5 

1.122 

1.120 

1.122 

1.123 

1.120 

1.120 

1.121 

1.120 

1.120 

1.120 

1.120 

N1-N5 

3.003 

4.165 

3.104 

2.905 

4.670 

4.102 

4.084 

4.095 

2.974 

4.139 

4.098 

N1-N2-N3 

156.9 

163.2 

157.6 

154.2 

114.5 

160.3 

165.1 

157.0 

155.3 

165.7 

157.8 

N2-N3-N4 

89.3 

72.0 

90.0 

93.6 

74.6 

78.9 

69.1 

81.9 

87.2 

73.4 

80.8 

N3-N4-N5 

101.2 

158.9 

105.7 

94.2 

170.9 

162.5 

169.5 

170.2 

18.9 

166.0 

166.9 

N4-N5-N1 

101.9 

154.4 

95.1 

109.6 

34.2 

21.0 

35.1 

26.5 

159.0 

27.2 

24.2 

N2-N1-N5 

88.8 

75.5 

87.4 

88.0 

86.6 

81.8 

80.0 

83.6 

77.5 

78.3 

83.0 

N2-N3-R 

129.6 

111.1 

136.3 

136.1 

113.1 

115.3 

127.1 

119.2 

123.3 

134.1 

116.4 

N4-N3-R 

134.9 

78.7 

129.2 

129.2 

38.5 

109.9 

73.9 

96.4 

147.4 

83.1 

100.0 

"  For  numbering  system,  see  Figure  1 . 

TABLE  11:  Monosubstituted  l-H,2-R-l,2-Pentazole  Decomposition  Product  Geometries  (Bond  Lengths  A— B  in  A  and  Angles 
A— B— C  in  Degrees)" 


H 

F 

ch3 

CN 

nh2 

OH 

CH2NO2 

n3 

nf2 

c2h3 

och3 

N2-R 

1.025 

1.358 

1.480 

1.336 

1.247 

1.346 

1.389 

1.132 

1.523 

1.425 

1.329 

N1-N2 

1.184 

1.235 

1.181 

1.190 

1.376 

1.206 

1.268 

1.323 

1.196 

1.183 

1.210 

N2-N3 

4.071 

5.719 

4.107 

3.961 

3.916 

4.066 

4.252 

4.071 

3.965 

4.078 

4.051 

N3-N4 

1.120 

1.120 

1.120 

1.120 

1.120 

1.120 

1.120 

1.139 

1.120 

1.120 

1.120 

N4-N5 

3.000 

2.797 

3.047 

2.946 

4.792 

3.071 

3.158 

1.296 

2.967 

3.044 

3.102 

N1-N5 

1.184 

1.168 

1.193 

1.188 

1.246 

1.195 

1.145 

1.431 

1.191 

1.197 

1.203 

N1-N2-N3 

83.2 

28.0 

82.2 

86.5 

75.6 

83.4 

115.4 

43.8 

86.6 

83.0 

85.4 

N2-N3-N4 

16.1 

15.4 

13.3 

18.6 

7.9 

12.8 

17.5 

23.4 

18.4 

12.7 

13.2 

N3-N4-N5 

156.1 

175.1 

152.4 

158.7 

32.9 

151.5 

163.5 

167.9 

158.4 

151.7 

149.6 

N4-N5-N1 

77.8 

122.4 

76.4 

78.3 

83.5 

74.8 

73.9 

106.8 

77.8 

76.0 

75.0 

N2-N1-N5 

157.6 

164.3 

158.5 

155.0 

114.9 

160.3 

169.5 

110.9 

155.3 

157.8 

157.6 

N1-N2-R 

129.2 

109.9 

135.6 

136.0 

113.3 

115.2 

115.4 

173.6 

123.2 

138.4 

116.4 

N3-N2-R 

140.3 

114.8 

135.2 

136.3 

37.7 

111.1 

142.9 

131.4 

148.3 

132.6 

105.1 

"  For  numbering  system,  see  Figure  1 . 


TABLE  12:  Activation  Energies  ( Ea )  for  the  Decomposition  of  Protonated  Pentazole  To  Form  the  Corresponding  Azidinium 
Cations  and  N2  (Including  Zero-Point  Energy  Correction,  in  kcal/mol)" 

1 ,3-Isomer  1 ,2-Isomer 

Ea  Ea\  Ea  2 


MP2 

CR-CCSD(T) 

MP2 

CR-CCSD(T) 

MP2 

CR-CCSD(T) 

H 

24.2 

28.2 

30.3 

30.9 

46.8 

46.1 

F 

22.6 

24.2 

19.3 

22.1 

36.7 

37.5 

ch3 

24.0 

28.8 

37.0 

35.6 

56.9 

53.9 

CN 

19.0 

22.4 

19.0 

22.3 

37.3 

39.3 

nh2 

18.7 

22.6 

30.5 

24.6 

50.5 

57.8 

OH 

19.9 

22.0 

23.6 

24.6 

33.5 

40.2 

CH2NO2 

22.6 

26.8 

28.6 

30.7 

17.9 

13.6 

n3 

16.0 

20.2 

23.4 

25.2 

14.1 

6.4 

nf2 

19.7 

23.3 

21.0 

23.4 

34.2 

39.3 

c2h3 

21.8 

26.8 

32.0 

31.1 

42.4 

50.8 

OCH3 

18.7 

21.9 

25.7 

26.8 

3.2 

16.0 

"Subscripts  1  and  2  label  the  first  and  second  activation  energies,  respectively,  for  the  1,2-isomer.  Ea2  is  calculated  relative  to  the  global  minimum 
energy  structures  shown  in  Figure  4. 


anions,  respectively.  The  differences  in  the  heats  of  reaction 
between  the  ionic  dimers  of  the  two  isomers  are  at  least  2.0, 
3.5,  and  6.5  kcal/mol,  respectively,  for  the  dinitramide,  nitrate, 
and  perchlorate  anions.  Although  the  parent  isomer  I  is  more 
stable  than  isomer  II,  the  anions  seem  to  prefer  to  bind  to  the 
latter. 

Conclusions 

Ionic  liquid  dimers  composed  of  the  five-membered  mono¬ 
substituted  pentazole  cation  and  the  dinitramide,  nitrate,  and 
perchlorate  anions  have  been  investigated.  Of  the  two  mono¬ 
substituted  pentazole  cation  isomers,  it  is  predicted  that  the  1,3- 
isomer  is  lower  in  energy  than  the  1,2-isomer.  The  relative 
activation  energies  and  decomposition  energies  for  the  decom¬ 
position  of  these  cations  to  N2  plus  the  corresponding  azidinium 
cation  can  be  related  to  the  electron-donating  character  of  the 


TABLE  13:  Reaction  Energies  Including  Zero-Point  Energy 
Correction  for  the  Decomposition  of  Protonated  Pentazole 
To  Form  the  Corresponding  Azidinium  Cations  and  N2 
(kcal/mol) 

1,3-Isomer  1,2-Isomer 

MP2  CR-CCSD(T)  MP2  CR-CCSD(T) 


H 

-18.3 

-4.0 

-38.3 

-23.8 

F 

-18.2 

-19.1 

-43.3 

-43.4 

ch3 

-16.7 

-3.0 

-35.1 

-21.0 

CN 

-21.8 

-7.1 

-42.5 

-26.5 

nh2 

-24.5 

-31.1 

-40.7 

-46.4 

OH 

-20.1 

-22.5 

-37.0 

-40.3 

CH2NO2 

-22.8 

-17.9 

-56.1 

-73.2 

n3 

-24.0 

-25.2 

-58.9 

-56.3 

nf2 

-20.0 

-3.5 

-38.0 

-22.1 

c2h3 

-21.0 

-17.9 

-35.4 

-21.8 

och3 

-19.9 

-22.2 

-35.0 

-38.1 
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(A) 


,  N 


(B) 


Figure  14.  Decomposition  pathways  for  the  monosubstituted  1,3- 
isomer  (A)  and  1,2-isomer  (B). 


.*%* 

0.0  kcal/mol  0.4  kcal/mol 

{0.0  kcal/mol)  (1.0  kcal/mol) 

Figure  15.  Minimum  energy  neutral  dimer  pairs  formed  by  the  pro- 
tonated  1,2-pentazole  (left)  and  1,3-pentazole  (right)  with  the  dinitra- 
mide  anion.  Relative  energies  (beneath  each  structure)  are  from  RHF/ 
6-311++G(d,p).  The  MP2/6-311++G(d,p)  relative  energies  at  RHF 
geometries  are  in  parenthesis.  Atom  colors  are  defined  in  Figure  3. 


TABLE  14:  Heats  of  Proton-Transfer  Reaction,  in  kcal/mol, 
from  the  Protonated  Pentazole  to  the  Anion  =  N(N02)2~, 
N03  ,  and  C104  a 


N(N02)2- 

NO3- 

CIO4- 

MP2 

1,3-Isomer 

-131.3 

-143.7 

-121.0 

CR-CCSD(T) 

-126.9 

-142.3 

-117.9 

MP2 

1,2-Isomer 

-150.7 

-163.0 

-140.3 

CR-CCSD(T) 

-145.8 

-161.3 

-136.9 

a  The  term  pentazolium  is  not  used  because  it  refers  to  the 
deprotonated  pentazole  (N5+). 


substituent.  An  electron- donating  substituent  tends  to  increase 
the  activation  energy  and  decrease  the  decomposition  energy. 
The  MCSCF  jt  orbital  analysis  suggests  that  the  electrons  in 
the  cationic  pentazole  ring  are  delocalized.  The  calculated 
enthalpies  of  formation  show  that  electron-withdrawing  groups, 
such  as  — CN  and  —  N3,  help  to  increase  the  heats  of  reaction. 
Several  dimer  structures  composed  of  the  H-substituted  penta¬ 
zole  and  either  of  the  anions  mentioned  above  are  proposed 
and  found  to  be  very  unstable.  The  low  barrier  connecting  the 


29.5 

33.7 

-110.4 


34.1 

-110.1 


Figure  16.  Ionic  pairs  in  the  parent  1,3-pentazole  and  dinitramide 
system  located  at  the  RHF/6-311++G(d,p)  level.  The  first  and  second 
energy  values  are  from  RHF/6-311++G(d,p)  and  MP2/6-311++G- 
(d,p)  at  RHF  geometries  in  kcal/mol.  The  third  energy  value  is  the 
MP2  heat  of  reaction  for  HN5H+  +  Y“  —  [HN5][HY]  (or  [HN5H+][Y-]), 
where  Y  =  N(N02)2,  NO3,  or  CIO4,  in  kcal/mol.  The  zero  of  energy 
for  Figures  16— 21  is  a  neutral  complex  presented  in  Figure  21.  Atom 
colors  are  defined  in  Figure  3. 


31.7  39.8 


Figure  17.  Transition  states  in  the  protonated  1,3-pentazole  and 
dinitramide  system,  connecting  ionic  pairs  to  neutral  complexes.  The 
energy  values  and  atom  colors  have  the  same  meaning  as  defined  in 
Figure  16.  Activation  energies  ( Ea )  from  RHF  and  MP2//RHF  are 
calculated  as  the  difference  between  the  transition  state  and  their 
corresponding  ionic  dimer  energies  in  kcal/mol. 

ionic  and  neutral  dimers  and  the  large  corresponding  exother- 
micity  suggests  that  a  proton  is  spontaneously  transferred  to 
the  anion.  So,  a  fundamental  step  for  the  decomposition  of 
protonated  pentazole  is  likely  to  be  deprotonation.  However, 
the  introduction  of  additional  ion  pairs  are  likely  to  stabilize 
the  charge  separation  in  the  ionic  species.  It  is  predicted  that 
the  1 -H,  3 -H-pentazolium— dinitramide  pair  with  a  barrier  for 
proton  transfer  of  4.7  kcal/mol  is  the  most  stable  dimer  among 
all  of  the  ion-pairs  studied.  The  ion-pair  comprised  of  the 
pentazolium  cation  and  perchlorate  anion  possesses  the  largest 
energy  content  and  would  have  been  a  desirable  high-energy 
ionic  liquid  but  the  2.7  kcal/mol  barrier  for  proton  transfer  is 
prohibitively  small  for  the  synthesis  and  detection  of  this  ionic 
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-137,1  -130.2 

Figure  18.  Neutral  dimers  in  the  protonated  1,3-pentazole  and 
dinitramide  system,  resulting  from  proton  transfers.  The  energy  values 
and  atom  colors  have  the  same  meaning  as  defined  in  Figure  16. 


0.0 

4.9 

10.3 

0.0 

4-9 

12.7 

-144.1 

-139.3 

-131.4 

Figure  21.  Neutral  dimers  in  the  parent  1,2-pentazole  and  dinitramide 
system,  resulting  from  proton  transfers.  The  energy  values  and  atom 

colors  have  the  same 

meaning  as  defined  in 

Figure  16.  The  lowest 

energy  dimer  found  is  the  left  structure  and  is  taken  as  the  energy  zero 

for  Figures  16—21. 
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Figure  19.  Ionic  pairs  in  the  parent  1,2-pentazole  and  dinitramide 
system.  The  energy  values  and  atom  colors  have  the  same  meaning  as 
defined  in  Figure  16. 


33.0  33.0  30.9 

25.4  263  26.8 


Figure  20.  Transition  states  in  the  protonated  1,2-pentazole  and 
dinitramide  system  connecting  ionic  pairs  to  neutral  complexes.  The 
energy  values  and  atom  colors  have  the  same  meaning  as  defined  in 
Figure  16.  Ea  values  are  defined  in  Figure  17. 


liquid.  The  l-H,3-H-pentazolium— dinitramide  pair  is  more  ideal 
with  a  relatively  larger  barrier  and  energy  content.  Among  all 
the  ion-pairs  studied,  it  is  calculated  that  the  pentazolium— nitrate 
pairs  provide  the  largest  energy  released  with  a  heat  of  reaction 
of  about  130  kcal/mol. 

The  CR— CCSD(T)  activation  energies  for  decomposition  of 
the  1,2-  and  1,3-  cations  are  in  the  range  of  22—31  and  20—28 
kcal/mol,  respectively.  Because  a  higher  barrier  corresponds  to 
a  more  stable  species,  the  most  desirable  substituents  from  this 
perspective  are  H,  CH3,  CH2NO2,  and  C2H3  for  the  both  isomers 
and  OCH3  as  well  for  the  1,2-isomer.  Because  all  of  the 
substituted  species  are  highly  exothermic,  the  stability  to 
decomposition  should  be  the  overriding  factor. 


ion  pair  t.s  neutral  pair 

Figure  22.  Structures  for  the  parent  1,3-pentazole  (upper)  and  1,2- 
pentazole  (lower)  with  nitrate:  ion  pair,  transition  state,  and  neutral 
pair.  The  energy  values  and  atom  colors  have  the  same  meaning  as 
defined  in  Figure  16.  Ea  values  are  defined  in  Figure  17. 
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Figure  23.  Structures  for  the  parent  1,3-pentazole  (upper)  and  1,2- 
pentazole  (lower)  with  perchlorate:  ion  pair,  transition  state,  and  neutral 
pair.  The  energy  values  have  the  same  meaning  as  defined  in  Figure 
16.  Atom  colors  are  defined  in  Figure  2.  Ea  values  are  defined  in 
Figure  17. 
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